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ABSTRACT 

We study the global stability of non-axisymmetric p-modes (also called inertial-acoustic 
modes) trapped in the inner-most regions of accretion discs around black holes. We 
show that the lowest-order (highest-frequency) p-modes, with frequencies lu = (0.5 — 
0.7)mQigcQ (where m = 1, 2,3, • • ■ is the azimuthal wave number, r^isco is the disc 
rotation frequency at the Inner-most Stable Circular Orbit, ISCO), can be overstable 
due to general relativistic effects, according to which the radial epicyclic frequency n is 
a non-monotonic function of radius near the black hole. The mode is trapped inside the 
corotation resonance radius Tc (where the wave pattern rotation speed uj/m equals the 
disc rotation rate fl) and carries a negative energy. The mode growth arises primarily 
from wave absorption at the corotation resonance, and the sign of the wave absorption 
depends on the gradient of the disc vortensity, ( = k^/ (20.11) (where S is the surface 
density). When the mode frequency oj is sufficiently high, such that dC/r > at rc, 
positive wave energy is absorbed at the corotation, leading to the growth of mode 
amplitude. The mode growth is further enhanced by wave transmission beyond the 
corotation barrier. We also study how the rapid radial inflow at the inner edge of the 
disc affects the mode trapping and growth. Our analysis of the behavior of the fluid 
perturbations in the transonic flow near the ISCO indicates that, while the inflow 
tends to damp the mode, the damping effect is sufficiently small under some condi- 
tions (e.g., when the disc density decreases rapidly with decreasing radius at the sonic 
point) so that net mode growth can still be achieved. We further clarify the role of 
the Rossby wave instability and show that it does not operate for black hole accretion 
discs with smooth- varying vortensity profiles. Overstable non-axisymmetric p-modes 
driven by the corotational instability provide a plausible explanation for the high- 
frequency (J> 100 Hz) quasi-periodic oscillations (HFQPOs) observed from a number 
of black-hole X-ray binaries in the very high state. The absence of HFQPOs in the soft 
(thermal) state may result from mode damping due to the radial infall at the ISCO. 
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1 INTRODUCTION 

1.1 Models of Quasi-Periodic X-Ray Oscillations: A Brief Review 

Rapid X-ray variabilities from Galactic compact binary systems have been studied for decades (e.g. van del Klis 2006). In recent 
years, our knowledge of quasi-periodic oscillations (QPOs) in black-hole X-ray binaries has greatly improved (see Remillard & 
McClintock 2006 for a review), thanks in large part to NASA's Rossi X-ray Timing Explorer (Swank 1999). The low- frequency 
QPOs (about 0.1-50 Hz) are common, observable when the systems are in the hard state and the steep power- law state (also 
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called "very high state"; see Done, Gierlinski & Kubota 2007), and they typically have high amplitudes and high coherence 
(Q > 10), and can vary in frequency on short timescales (minutes). However, it is the weaker, transient High- Frequency QPOs 
(HFQPOs, 40-450 Hz) that have attracted more attention, since their frequencies do not vary significantly in response to 
sizable (factors of 3-4) luminosity changes and are comparable to the orbital frequencies at the Innermost Stable Circular 
Orbit (ISCO) of black holes with mass M ~ lOM©. As such, HFQPOs potentially provide a probe to study the effects of 
strong gravity. HFQPOs are usually observed in the very high state of the X-ray binaries, and have low amplitudes (0.5 — 2% 
rms at 2-60 keV) and low coherence (Q ~ 2 — 10). Out of the seven black-hole binaries from which HFQPOs have been 
reported, four show pairs of QPOs (first discovered in GRO J1655-40; Strohmayer 2001) with frequency ratios close to 2 : 3 
(300 and 450 Hz in GRO J1655-40, 184 and 276 Hz in XTE ,11550-564, 113 and 168 Hz in GRS 1915+105, 165 and 240 Hz 
in H1743-322; note that GRS 1915-1-105 also has a second pair of QPOs with / = 41 and 67 Hz). 

It is worth noting that QPO (with period of ~ 1 hour) in X-ray emission has recently been detected in the active galaxy 
RE ,11034+396 (Gierhnski et al. 2008). This could be the "supermassive" analog of the HFQPOs detected in black-hole X-ray 
binaries. 

Despite the observational progress, the origin of the HFQPOs remain unclear. A number of possibilities/models have 

been suggested or studied to various degrees of sophistication. We comment on some of these below: 

- Stella, Vietri & Morsink (1999) and others (see Schnittman & Bertschinger 2004, Schnittman 2005) suggested that 
orbiting hot spots (blobs) in the disc oscillating with epicyclic frequencies may provide variability in the X-ray emission. 
However the radial positions of such blobs are free parameters, which must be tuned to match the observed QPO frequencies, 
and it is also not clear that the blobs can survive the differential rotation of the disc. 

- Abramowicz & Kluzniak (2001) suggested that HFQPOs involve certain nonlinear resonant phenomenon in the disc 
(e.g., coupling between the radial and vertical epicyclic oscillations of the disc fluid element; Kluzniak & Abramowicz 2002). 
This was motivated by the observed stability of the QPO frequencies and the commensurate frequency ratio. However, so far 
analysis has been done based only on toy models involving coupled harmonic oscillators (e.g. Rebusco 2004; Horak & Karas 
2006) and no fluid dynamical model producing these resonances has been developed (sec Abramowicz et al 2007 and Rebusco 
2008 for recent reviews). Petri (2008) considered the resonant oscillation of a test mass in the presence of a spiral density 
wave, but the origin of the wave is unclear. 

- Acoustic oscillation modes in pressure-supported accretion tori have been suggested as a possible source of the observed 
QPOs (RezzoUa et al. 2003; Lee, Abramowicz & Kluziniak 2004; see also Blaes, Arras & Fragile 2006; Schnittman & RezzoUa 
2006, Blaes et al 2007, Sramkova et al 2007). In this model, the commensurate mode frequencies arise from matching the 
radial wavelength to the size of the torus. Note that the QPO frequencies are determined mainly by the radial boundaries 
of the torus, which must bo tuucd to match the observed QPO frequencies. It is also not clear that the accretion flow in the 
very high state (in which HFQPOs axe observed) is well represented by such a torus (e.g. Done et al 2007). 

- Li & Narayan (2004) considered the dynamics of the interface between the accretion disc and the magnetosphere of a 
central compact object (see also Lovelace & Romanova 2007) . The interface is generally Rayleigh-Tayor unstable and may also 
be Kelvin-Helmholtz unstable. While such an interface is clearly relevant to accreting magnetic neutron stars, Li & Naxayan 
suggested that it may also be relevant to accreting black holes and that the strongly unstable interface modes may give rise 
to QPOs with commensurate frequencies. 

- Perhaps the theoretically most appealing is the relativistic diskoseismic oscillation model, according to which general 
relativistic (GR) effects produce trapped oscillation modes in the inner region of the disc (Kato & Fukue 1980; Okazaki et 
al. 1987; Nowak & Wagoner 1991; see Wagoner 1999; Kato 2001 for reviews; see also Tassev & Bertschinger 2007 for the 
kinematic description of some of these wave modes) . A large majority of previous studies have focused on disc g- modes (also 
called inertial modes or inertial-gravity modes, whose wavefunctions - such as the pressure perturbation, contain nodes in the 
vertical direction), because the trapping of the g-mode does not require a reflective inner/outer disc boundary. Kato (2003a) 
and Li, Goodman & Narayan (2003) showed that the g-mode that contains a corotation resonance (where the wave patten 
frequency equals the rotation rate of the background flow) in the wave zone is heavily damped. Thus the only nonaxisymmetric 
(m 0) g-modes of interest are those trapped around the maximum of H + K/m (where O is the rotational frequency, k is 
the radial epicyclic frequency and in is the azimutlial mode number; see Fig. 1 below). Unfortunately, the frequencies of such 
modes, uj ~ mnisco, are too high (by a factor of 2-3) compared to the observed values, given the measured mass and the 
estimated spin parameter of the black hole (Silbcrglcit & Wagoner 2007; see also Tassev & Bertschinger 2007) . Axisymmetric 
g-modes (m = 0) may still be viable in the respect, and recent studies showed that they can be resonantly excited by global 
disc deformations through nonlinear effects (Kato 2003a,2008; Ferreira & Ogilvie 2008). Numerical simulations (Arras, Blaes 
& Turner 2006; Reynolds & Miller 2008), however, indicated that while axisymmetric g-mode oscillations are present in the 
hydrodynamic disc with no magnetic field, they disappear in the magnetic disc where MHD turbulence develops. Also, Fu & 
Lai (2008) carried out an analytic study of the effect of magnetic fields on diskoseismic modes and showed that even a weak 
(sub-thermal) magnetic field can "destroy" the self-trapping zone of disc g-modos, and this may (at least partly) explain the 
disappearance of the g-modes in the MHD simulations. 

- Tagger and collaborators (Tagger & Pellat 1999; Varniere & Tagger 2002; Tagger & Varniere 2006; see Tagger 2006 
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for a review) developed the theory of accretion-ejection instability for discs threaded by strong (of order or stronger than 
equipartition) , large-scale poloidal magnetic fields. They showed that such magnetic field provides a strong coupling between 
spiral density waves and Rossby waves at the corotation, leading to the growth of the waves and energy ejection to disc corona. 
Tagger & Varniere (2006) suggested that normal modes trapped in the inner region of the disc become strongly unstable by a 
combination of accretion-ejection instability and an MHD form of the Rossby wave instability (see Lovelace et al. 1999; Li et 
al. 2000; see section 6 below). The Tagger model has the appealing feature that the instability leads to energy ejection into 
the disc corona, and thus explains why HFQPOs manifest mainly as the variations of the nonthermal (power-law) radiation 
from the systems. 



1.2 This Paper 

In this paper we study the global corotational instability of nonaxisymmetric p-modes (also called inertial-acoustic modes) 
trapped in the inner-most region of the accretion disc around a black hole. The p-modes do not have vertical structure (i.e., 
the wavefunctions have no node in the vertical direction). We focus on these modes because their basic wave properties (e.g. 
propagation diagram) are not affected qualitatively by disc magnetic fields (Fu & Lai 2008) and they are probably robust 
under hydromagnetic effects and disc turbulence (see Reynolds & Miller 2008). 

The corotational instability of p-modes studied in this paper relies on the well-known GR effect of test-mass orbit around 
a black hole: Near the black hole, the radial cpicyclic frequency k reaches a maximum (at r = 8GM/c^ for a Schwarzschild 
black hole) and goes to zero at the ISCO (risco = QGM/c?). This causes non-monotonic behavior in the fluid vortensity, 
C = re^/(2E0) (assuming the surface density E is relatively smooth), such that dC,/dr > for r < rpeak (where Tpeak is the 
radius where <^ peaks) and dC^/dr < for r > Tpoak- The vortensity gradient dC,/dr plays an important role in wave absorption 
at the corotation resonance (Tsang & Lai 2008a; see also Goldreich & Tremaine 1979 for corotational wave absorption due to 
external forcing) . We show that the p-modes with frequencies such that the corotation radii lie inside the vortensity peak can 
grow in amplitude due to absorption at the corotation resonance, and that the ovcrstability can be achieved for several modes 
with frequencies closely commensurate with the azimuthal wavenumber m. Tagger & Varniere (2006) have studied similar 
modes in discs threaded by strong magnetic fields, but in our analysis magnetic fields play no role. 

The trapping of the p-modes requires the existence of a (partially) reflecting boundary at the disc inner edge, close to 
the ISCO. One may suspect that the rapid radial inflow at the ISCO will diminish any potential instabilities in the inner 
accretion disc (see Blaes 1987 for the case of thick accreting tori). Our analysis of the wave perturbations in the transonic 
accretion flow (sec section 5) suggests that waves are partially reflected at the sonic point, and global overstable p-modes may 
still be produced under certain conditions (e.g., when the surface density of the flow varies on sufficiently small length scale 
around the sonic point). Even better mode trapping (and therefore larger mode growth) may be achieved when the system is 
an accretion state such that the inner disc edge does not behave as a zero-torque boundary (see section 7 for discussion and 
references). 

Our paper is organized as follows. After summarizing the basic fluid equations for our problem in section 2, we give 
a physical discussion of the origin of the corotational instability of disc p-modos in section 3. We present in section 4 our 
calculations of the growing p-modes with simple reflective inner disc boundary conditions. Section 5 contains our analysis of 
the effect of the transonic radial inflow at the ISCO on the p-mode growth rate. In section 6, we discuss the role of the Rossby 
wave instability and show that it is not effective in typical accretion discs under consideration. In section 7 we discuss the 
application of our results to HFQPOs in black hole X-ray binaries. 



2 SETUP AND BASIC EQUATIONS 



Wo consider a geometrically thin disc and adopt cylindrical coordinate system {r,(f>,z). The flow is assumed to bo barotropic, 
so that the vertically integrated pressure, P = J p dz, depends only on the surface density, E = J p dz. We use the pseudo- 
Newtonian potential of Paczynski & Wiita (1980) 

* = -^, (1) 

r -rs 

with rs = 2GM/(? the Schwarzschild radius. The free-particle (Keplerian) orbital and radial cpicyclic (angular) frequencies 
are 



r dr 



1/2 



GM 



rs 



—(r fix 

r dr 



1/2 



'irs 



■ rs 



(2) 



The function k peaks at r = (2 -|- ^/S)rs and declines to zero at nsco = 3rs (while for a Schwarzschild black hole in GR, 
K peaks at r = 4rs). The unperturbed flow has velocity uo = (ur,rQ,,0). Since pressure is negligible for thin discs, we have 
O ~ Qk. 
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Neglecting the self-gravity of the disc we have the linear perturbation equations: 

^5E + V- (E5u + uo(5E) =0, (3) 

^(5u+ (uo • V)(5u+ (<5u- V)uo = -V(5/i, (4) 

where (5E, 5\i and 5h = 5P/E are the (Eulerian) perturbations of surface density, velocity and enthalpy, respectively. For 
barotropic flow, 5h and 5E are related by 

^E 

&h = cl — , (5) 

where is the sound speed, with = dP/d'E. We assume all perturbed quantities to be of the form e*""^"'"*, where m is a 
positive integer, and uj is the wave (angular) frequency. The perturbation equations then become 

-iCb^Sh+'^CErSur) + —T,Su4, + --^(rurSE) = 0, (6) 
ci r or r r or 

d d 

— iCbSUr — 2Q,5us + — {UrSUr) = ——5h, (7) 

or or 

—lUjdUA, + ——dUr ^ — {rbuA,) = oh, (8) 

2\l r or r 

where 

Cii = Lu — mO,, (9) 

is the wave frequency in the frame corotating with the unperturbed fluid. 

Except very near the inner edge of the disc, rin — risco, the unperturbed radial velocity is small, \ur\ <^ rO,. In our 
calculations of the disc modes, we will neglect Ur and set the last terms on the left-hand sides of equations ©-(O to zero 
(However, Ur plays an important role in determining the inner boundary condition of the fluid perturbations at rin; see section 
5). Eliminating the velocity perturbations in favor of the enthalpy, we obtain our master equation 

d^ d f, D\ d 2mn f d , ^2E^ D] 

oh — 0, (10) 



dr^ dr \ rE J dr ruj \ dr D 



where 

D = K-Co^. (11) 
For concreteness we assume the surface density to have a the power-law form 

Eocr-f, (12) 

where p is the density index. 

The above equations adequately describe disc p-modes (also called inertial-acoustic modes), which do not have vertical 
structure. Other disc modes (g-modes and c-modes) involve the vertical degree of freedom [see Kato 2001 for a review; also 
see Fig. 1 of Fu & Lai (2008) for a quick summary], and their stability properties are studied by Kato (2QQ3a), Li et al. (2003) 
and Tsang & Lai (2008b). 

To determine the global modes of the disc, appropriate boundary conditions must be specified. These are discussed in 
sections 4 and 5. 



3 P-MODES AND THEIR GROWTH DUE TO COROTATION RESONANCE: A PHYSICAL 
DISCUSSION 

A WKB analysis of the wave equation (|10p yields the dispersion relation for the local plane wave 5h oc exp [i J^k{s)ds\ : 

,2 D 2mQ, / d nE\ 

^ +72+^h7zln^h^0. (13) 



ruj \dr D 

Far from the singularity {uj — 0) at the corotation radius rc, this reduces to the well-known dispersion relation of spiral density 
wave with no self-gravity (e.g., Shu 1992), fe^ ~ —D/Cg, or 

Cj'^ ^ k'^ + k^cl . (14) 

Density waves (p-modes) can propagate inside the inner Lindblad resonance radius riL (defined hj u = — k), and outside the 
outer Lindblad resonance radius roL (defined by iIj = k), i.e., in the region where uo/m < f2 — «:/m and u/m > Q. -\- K/m, 
respectively (see Fig. 1). Between riL and roL, waves are evanescent except that a very narrow Rossby wave zone exists around 
the corotation radius. Indeed, in the vicinity of a) = 0, equation (|13p reduces to 

2mQ, f d\nC,\ , , 

(15) 
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radius r 

Figure 1. Wave propagation diagram for non-axisymmctric p-modcs in thin accretion discs around black holes. In the upper panel, 
the three solid curves depict the disc rotation profile and H ± K/m (where k is the radia epicyclic frequency); note that the three 
curves join each other at the disc inner radius rin = nsco (the inner-most circulaj: orbit) since /t(r'isco) = 0- The wavy lines (of height 
oj/m) indicate the propagation zones for inertial-acoustic waves. Disc p-modes are trapped between r-jn and the inner Lindblad resonance 
radius (where uj/m = CI — K/m), but can tunnel through the corotation barrier. The lower panel depicts the disc vortensity profile, 
= K^/(2r2S), which has a maximum at the radius rpoak (^s long as the surface density S docs not vary too strongly with r). P-modcs 
with Lu/m > Hpcaij = n(rpoak) (the upper wavy line) arc ovcrstablc due to wave absorption at the corotation resonance radius rc (where 
Lu/m = Q) since {dC/dr)r^ > 0. P-modcs with cj/m < fipoak (the lower wavy line) tend to be damped by wave absorption at Tc since 
{d(^/dr)r^ < 0. Note that a narrow Rossby wave zone (labeled by thick horizontal bars) exists just outside or inside the corotation radius 
— the location of this Rossby zone determines the sign of the corotational wave absorption. 



where 



c = 



2nE 



(16) 



is the vortensity of the (unperturbed) flow. For {dC,/dr)r^ > 0, the Rossby wave zone lies between rc and Tc + Ar^, where 
Arfl = {2cs/k)\v\ and p^^'''' kdr = -Kv, with the number of wavelengths in the Rossby zone given by (Tsang & Lai 2008a) 



Cs dlnC 
QK dr 



dr 



r 



(17) 



where q = — (dln^l/dlnr)^ , and the second equality assumes E oc r"*". For {d(/dr)rc < 0, the Rossby wave zone lies inside 
Fc, between Tc — |Arij| and rc (see Fig. 1). Note that since v ~ Cs/{Kr) ~ H/r <^ 1, no standing Rossby wave can exist in the 
Rossby zone (see also section 6). 

Assuming that there exists a reflecting boundary at the inner disc radius rin ^ risco (see sections 4.3 and 5), normal 
modes can be produced, with the waves partially trapped between rin and riL - these are the p-modes that we will focus on 
in this paper. The mode eigen-frequency u; = u>r + iiOi is generally complex, with the real part ujr determined approximately 
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by the Sommerfeld "quantization" condition 



''''' 1 / 

— V iDr — K^dr = mx + ip, (18) 



where ujr = uOr — mO,, n is an integer and ip (of order unity) is a phase factor depending on the details of the (inner and outer) 
boundary conditions. The overstability of the p-mode is directly related to the reflectivity of the corotation barrier between 
riL and roL- In the WKB approximation, the imaginary part of the mode frequency is given by (Tsang & Lai 2008a; see also 
Narayan et al. 1987, who considered shearing- sheet model) 

where TZ is the reflectivity (see below). Thus the mode becomes overstable [uoi > 0) for \'JZ\ > 1 (termed "super-reflection") 
and stable [uji < 0) for \n\ < 1. 

Super-reflection in fluid discs arises because the waves inside the corotation radius and those outside carry energy or 
angular momentum of different signs: Since the wave inside rc has pattern speed uir/m less than the fluid rotation rate 
f2(r), it carries negative energy; outside rc, we have ujr/m > ri{r), the wave carries positive energy. Consider an incident 
wave Sh (x exp(— i kdr), carrying energy of the amount (—1), propagating from small radii toward the corotation barrier 
0. The wave reflected at riL takes the form Sh cx 7?.exp(i J' kdr), and the transmitted wave in the region r > roL is 
Sh oc Texp(j kdr). Because of the corotation singularity, the wave energy can also be transferred to the background flow 
and dissipated at the corotation radius. Energy conservation then gives —1 — ( — l)|7?.p + \T\^ + T)c, or 

|7^|^ = l + |T|^+Dc, (20) 

where Dc is the wave energy dissipated at the corotation. 

Tsang & Lai (2008a) derived the analytical expressions (in the WKB approximation) for T, TZ and "Dc- Two effects 
determine the reflectivity, (i) The transmitted wave (corresponding to the |Tp term) always carries away positive energy and 
thus increases |7?.|^. (ii) Wave absorption at the corotation can have both signs, depending on v: For v > Q, the Rossby wave 
zone lies outside rc, positive wave energy is dissipated and we have T>c > 0; for i/ < 0, the Rossby zone lies inside Vc and we 
have Oc < (see Fig. 1). Tsang & Lai (2008a) showed explicitly that under most conditions, \Vc\ 2> jTp (except when u ~ 0, 
for which T>c — 0). In the limit of ^ 1, we have 

|T|' ~ exp(-2eii), Oc ~ 27ri/exp(-2eiia), (21) 



where 



e„.p^^, e„..r^^^ (22) 

[these expressions are valid for On, Ona 2> 1; see Tsang & Lai (2008a) for more general expressions]. Thus super-reflectivity 
{\TZ\^ > 1) and growing modes (uji > 0) are achieved when 



1 /■'■OL _ 

> J'crit = --r- exp(-20iib), with 6iib = 0ii - Ona = / -■ (23) 

Note that typically ji^critl <C 1; if is neglected compared to Vc, then fcrit = 0. 

As mentioned before, since k, is non-monotonic near the black hole, the vortensity ( is also non-monotonic, attaining a 
peak value at r = rpeak before dropping to zero at the ISCO. Therefore, p-modes with frequencies such that the corotation 
radius lies inside rpoak are expected to be overstable by the corotational instability discussed above. In other words, when 
uJr/m > fipoak = ^^(rpoak), the corotatiou resonance acts to grow the mode. Note that fipcak depends on the surface density 
profile as well as the spacetime curvature around the black hole (see Fig. 2). On the other hand, when cur/m < fipeak 
(i' < 0), the corotational wave absorption acts to damp the mode. However, when ujr/m is only slightly smaller than ilpeak 
(z^crit < V < 0) mode growth can still be obtained due to wave leakage beyond the outer Lindblad resonance, though the 
growth rate will be small (see section 4.4 for examples). 



4 CALCULATIONS OF TRAPPED, OVERSTABLE P-MODES 

To determine the eigenvalues ujr and uji of the trapped modes, we solve equations ©-([H]) (with u,- = 0) or equation (flO 
subjected to appropriate boundary conditions at rin and rout- 



^ Note that since the group velocity of the wave has opposite sign as the phase velocity for r < rjL, the wave of the form exp(— j k dr) 
(with fc > 0) is outward propagating. 
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T 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




-1 1 

density index p 



Figure 2. Critical mode frequency for corotational instability as a function of the disc surface density index p (with S oc r~P). Wave 
absorption at the corotation resonance acts to grow the mode only if the corotation occurs in the region of positive vortensity gradient, 
i.e., if the mode pattern frequency Lj/m > Opcak (see Fig. 1). P-mode trapping also requires Lo/m < S^igco = ^(nsco)- Note that weak 
mode growth can still occur when uj/m is slightly below Op^ak due to wave leakage beyond the outer Lindblad resonance radius. See 
text for detail. 



4.1 "Landau" Integration Contour 

When solving eigenvalue problem using the standard method (e.g. the shooting method as described in Press et al 1998), 
we encountered a conundrum: For v > 0, we could find both a growing mode and a decaying mode, with almost the same 
ujr but opposite uji. This appears to contradict our discussion in section 3. This conundrum arises because our numerical 
integration is confined to the real r axis. However, analogous to Landau's analysis of wave damping in a plasma (e.g., Lifshitz 
& Pitaevskii 1981), care must be taken in defining appropriate contour of integration across the corotation resonance. Indeed, 
at corotation, equation (|10p contains a singular term, proportional to 



where Rc = rc — ircCUi / (qur) is the complex pole, rc is determined by uir ~ mQ{rc) and q — ~{dlnQ/d\nr)c > 0. 

As discussed in Lin (1955) in the context of hydrodynamical shear flows, to obtain physically relevant solutions of the 
fluid system, it is necessary that the integration contour lies above the pole. This is the Landau contour. In essence, only 
by adopting such a Landau contour can one obtain the correct wave absorption (dissipation) at the corotation. For growing 
modes {cui > 0), Im(i?c) < 0, our numerical integration along the real r axis constitutes the correct Landau contour. On the 
other hand, for decaying modes {uji < 0), the real r axis is not the correct Landau contour as Im(i?c) > 0. Instead, to obtain 
physical solutions for these decaying modes, the integration contour must be deformed so that Rc lies below it (see Fig. 3). 
As we are primarily interested in over-stable modes in this paper, it is adequate to integrate along the real r axis in our 
calculation. 
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Re(r) 

• Rc(cA>0) 


Im(i-) 

(5 
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^ Re(r) 



Figure 3. "Landau" contour for integration across the corotation resonance. To calculate the growing mode {u)i > 0), it is adequate to 
integrate the fluid perturbation equations along the real r axis (upper panel). To obtain the physical solution for the shrinking mode 
(u}i < 0), the integration contour must deformed so that the corotational pole Rc lies below the contour. 



4.2 Outer Boundary Condition 

As we are interested in self-excited modes in the inner region of the disc, we adopt the radiative outer boundary condition. 
Specifically, far from the outer Lindblad resonance (r > tol) we demand that only an outgoing wave exists: 

ShoaAexpfi ( k dr\ , with ^ = f^^) ' (25) 



where k = yJ-D/ci (see Tsang & Lai 2008a). This gives the boundary condition at sorriG /"out ^ ^OL • 

<5/i'(rout) = 5ft(r-out) (ifc + i^") . (26) 
In practice, we find that rout ~ 2roL would yield sufficiently accurate results. 



4.3 Inner Boundary Conditions 

To obtain global trapped modes, at least partial wave reflection must occur at rin. To focus on the effect of corotational 
instability discussed in section 3, in this section we consider two simple inner boundary conditions. We defer our analysis of 
the effect of radial inflow on the p-modes to section 5. 

(i) At the ISCO, the flow plunges into the black hole, we expect a sudden decrease in the surface density of the disc. 
Thus, it is reasonable to consider the free surface boundary condition, i.e., the Lagrangian pressure perturbation AP = 0. 
Using 5ur — —iili^r (where is the radial Lagrangian displacement), and AP = 5P + S^rdP/dr, we have 

where we have used dP/dr = [dP / dT,)(dTi/ dr) — —pclTi/r for barotropic, power-law discs (E oc r"''). 

(ii) We assume that the radial velocity perturbation vanishes at the inner boundary, i.e., 5ur — 0. This was adopted by 
Tagger & Varniere (2006) in their calculations of overstable global modes due to accretion-ejection instability. 

Both of these boundary conditions correspond to zero loss of wave energy at the inner boundary: If a wave from large 
radii impinges toward rin, the reflected wave will have the same amplitude. However, the phase shifts due to reflection differ 
in the two cases, and the resulting mode frequencies uj,. are different. Since the corotational wave amplification depends on 
uur, the mode growth rate uji will also be different. 
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Figure 4. Example wavefunctions for disc p-modes. The upper and middle panels show 5h and i&Ur (the solid lines for the real part 
and dashed lines for the imaginary part), the lower panels show the angular momentum flux, all in arbitrary units [with 5/i(rout) = 1]. 
The radius r is in units of GMji? . The disc sound speed is Cs = 0.1rf2, and the m = 2 modes are obtained using the inner boundary 
condition AP(r'igco) = 0- The left panels show the p-mode for the disc with a density profile S oc r~^, with the eigenvalues ujr = 
0.467mf2isco 1 ^il^v = 0.0029 [where f2isco = f^(^isco)]; the right panels show the mode for the disc with constant surface density 
profile (p = 0), with eigenvalues ojr = 0.464mf!isco i ^il^r = 0.00073. Note that the model shown on the left panels has Vc < Tpeak (the 
radius of peak vortensity) and thus F{rc+) < F{rc—), while the model shown on the right panels has rc > ?'peak ^^<i F{rc+) > F{rc—). 
In both models there is a positive flux for r > r^. due to the outward propagating wave. The inserts on the left panels show the blowups 
of the real wavefunctions near the corotation radius. 




density index p density index p 

Figure 5. The real and imaginary frequencies of disc p-modes (with azimuthal wave numbers m = 1,2,3) as a function of the surface 
density index p (where S oc r~P). The modes are calculated assuming the inner boundary condition AP(rjgco) = 0. The left panels are 
for discs with Cs = O.lrQ and the right panels for Cs = 0.2rr2. The dotted lines denote the lower bound ujr/m = Hpeak for which the 
corotational wave absorption acts to enhance mode growth. 
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Figure 6. The real and imaginary frequencies of disc p-modes (with m = 2,3) as a function of the normalized sound speed Cs/irQ). 
The disc is assumed to have a constant density profile (p = 0), and the inner boundary condition is AP(rjsco) = 0. The bottom panel 
shows the v parameters for the modes. 



4.4 Numerical Results 

We solve for the complex eigen-frequency = ojr + iuJi using the shooting method, with a fifth-order Runge-Kutta integrator 
(Press et al. 1992). As discussed in section 4.1 we only calculate the growing modes {uji > 0). We consider disc models with 
different surface density profile (characterized by the index p), sound speed Cs, and inner boundary conditions. For a given set 
of disc parameters and the azimuthal mode wavenumber m, the lowest order (highest frequency) mode has the best chance of 
being overstable. This is easily understood from our discussion in section 3 (see Fig. 1): a low- frequency wave has to penetrate 
a wider evanescent barrier for the corotational amplifier to be effective, and when lot < mfipoak the corotation resonance acts 
to damp the mode. For most disc models we have considered, the lowest order mode (of a given m) is the only mode that has 
ijJi > 0. 

Our numerical results are presented in Figures |4}{7l Figure 4 gives two examples of the eigenfunctions of overstable 
trapped p-modes, obtained with the inner boundary condition AP = 0. In addition to Sh and Sur, we also plot the angular 
momentum fiux carried by the wave across the disc (e.g., Goldreich & Tremaine 1979; Zhang & Lai 2006) 



F{r) = Trr^ERe {SurSi 



■KmrYl 



Im Sh 



dSh* 



(28) 



D \ dr 

where the second equality follows from equations (O-® (with Ur — 0). We see from Fig. 4 that outside the corotation radius 
(re), F is nearly constant since only the outgoing wave exists in this region and the wave action is conserved in the limit of 
uji <C LJr- Inside rc, the interference between the ingoing and outgoing waves gives rise to the variation of F. At risco, F 
approaches zero since no wave action is lost through the disc inner boundary when AP = 00 Figure 4 also shows a flux jump 
across the corotation. In the limit of uji <^ ujr, dSh/dr, Sur and are discontinuous across Vc (although Sh is continuous), 
giving rise to the flux discontinuity (see Tsang & Lai 2008a): 



f (re+) - F{r,-) 



sD 



\Sh\- 



(29) 



This discontinuity signifies the corotational wave absorption, the sign of which depends on oc {d(/dr)r^, as discussed in 
section 3. Thus, the model shown on the left panels of Fig. 4 has rc < r-peak and f > 0, and the mode growth is primarily 
driven by wave absorption at the corotation. The model shown on the right panels of Fig. 4, on the other hand, has > rpcak 
and 1/ < 0, thus the corotational wave absorption acts to damp the mode, and the overall growth of the mode is due to the 
outgoing wave beyond rc, and the growth rate is much smaller than the model shown on the left panels. 



^ Note that equation II28I I is the time averaged flux and is defined for waves with real u). Using equation II27I I and DSur = iuidSh/dr — 
(2imQ,/r)5h (obtained from eqs. [7]-[8]) it is easy to show that P(rin) = exactly for real uj. 
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Figure 7. The real and imaginary frequencies of disc p- modes (with azimuthal wave numbers m = 1,2, 3) as a function of the surface 
density index p (where S oc r~P). The modes are calculated assuming the inner boundary condition 5'!ir(''isco) = 0. The left panels 
are for discs with Cs = O.lrf! and the right panels for Ca = 0.2rf!. The dotted lines denote the lower bound uj/m = ^pcak for which the 
corotational wave absorption acts to enhance mode growth. 



Figures 5-6 show the frequencies of the fundamental (no node/highest frequency) growing p-modes (with m — 1,2,3) 
for various disc parameters, again obtained with the inner boundary condition AP = 0. We consider p in the range between 
— 1.5 and 1.5, and Cs up to O.Srfi. For a given sound speed, the real mode frequency ujr depends very weakly on p, but the 
growth rate uji increases with p (see Fig. 5) since a larger value of p leads to a larger v and enhanced wave absorption at the 
corotation [see equation p7[) ]. In general, as the sound speed increases, the effective wavelength of the mode increases, and 
LJr decreases in order "fit in" the trapping zone between rin and riL (see Fig. 6). The mode growth rates uji depends on Cs in 
a non-monotonic way because of two competing effects: As Cs increases, less attenuation occurs in the evanescent zone, and 
more wave energy can be absorbed at the corotation and propagate to the outer edge of the disc; these tend to increase uji. 
On the other hand, increasing Cs also leads to smaller cur, which shifts the corotation resonance to a larger radius and leads 
to decreasing v and uJi. 

Note that the growing modes shown in Fig. 5 extend below the ujr/m = S7poak boundary due to the propagation of waves 
beyond the corotation radius, as discussed in section 3 [see eqs. (|23|l ]. Such modes (with z/crit < < 0) grow significantly 
slower than the modes with u > as the flux is attenuated by the entire barrier between riL and roL- 

For comparison. Figure 7 shows the disc mode frequencies and growth rates when the inner boundary condition Su,. — 
is adopted. The different boundary condition leads to a different phase shift ip and higher mode frequency, but the results are 
similar to those illustrated in Fig. 5. In particular, as p increases, uJr remains approximately constant while cui increases. 



5 EFFECT OF RADIAL INFLOW ON THE P-MODE GROWTH RATE 

Our mode calculations presented in section 4 neglect the radial velocity of the accretion flow and assume a loss-less inner disc 
boundary condition (either AP = or 5ur = at risco)- In real discs, the radial inflow velocity Ur is not negligible as r 
approaches risco, and the flow goes through a transonic point (where Ur = —Ca) at a radius very close to risco- We expect 
that part of the fluid perturbations may be advected into the black hole and the inner disc boundary will not be completely 
loss-less. Here we study the effect of the transonic flow on the p-mode growth rate. 

We note that just as the accretion disc is not laminar but turbulent, the accretion flow around risco is complicated. 
General relativistic MHD simulations in 3D are only beginning to shed light on the property of the black hole accretion flow 
(e.g., Beckwith, Hawley & Krolik 2008; Shafee et al. 2008; Noble, Krolik & Hawley 2008), and many uncertainties remain 
unresolved. Here, to make analytic progress, we adopt a simple viscous transonic flow model, which qualitatively describes 
the inner accretion flow of the black hole as long as the flow remains geometrically thin (see Afshordi & Paczynski 2003). 
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5.1 Boundary Condition at the Sonic Point 

We rewrite equations ©-(IT]) as 



—rrSh' + 5u'^ = 
c 



Sh+ — Sur - — Su4, = Ai, (30) 

Ur r 



5h' + Ur Su'r — [iio — u[.) SUr + 2f2 Su^ = ^2, (31) 
where ' stands for d/dr and we have used rSur = constant for the background flow. Solving for Sh' and Su'^ we have 

5h' = :^-^, (32) 
1 - ui/cj 

(33) 

Clearly, in order for the perturbation to be regular at the sonic point Vs, where Ur = — Cs, we require 

A2 + CsAi =0 at r = rs. (34) 
For definiteness, we characterize the variations of E and at the sonic point — rjsco by the two length scales: 

=± =1. 

From rTjUr = constant, we also have itj, — Cs{r~^ + L^^) at r = r^. Then equation (|34|l becomes 



(35) 



icD 2 



Sh + 



iCu — 2cg 



1 1 

r Le 



r 



+ ( 2$! — 1 (5?i0 = at r = rs . (36) 



This is the boundary condition for the fluid perturbations at the sonic point. 
5.2 Properties of the Transonic Flow 

Before exploring the effect the radial inflow on the disc modes, we flrst estimate the length scale for the surface density 
variation, Ls, using the viscous slim disc model (e.g., Muchotrzeb & Paczynski 1982; Matsumoto et al. 1984; Abramowicz et 
al. 1988) 

The basic steady-state slim disc equations are 

M = -2nrT,Ur, (37) 
S' 

UrUj. = —Cs — + {fl ~ (38) 

2-1 

Mlo = Ml + 27ri/visr^En', (39) 

where I = r^Q, is the speciflc angular momentum of the flow, Q is the actual rotation rate, Qk is given by equation ((2)|, z^vis 
is the kinetic viscosity, and lo is the eigenvalue that must be solved so that flow pass through the sonic point smoothly. We 
shall use the a-disc model, so that i^vis ~ aHcs, with H ~ Cs/i^K- 

To estimate Le, we assume l{r) ~ r^QK{r) for r ^ nsco and lo — lK{risco)- Equation (|39|) gives 

uril - lo/l) ^ -3z/vi./(2r-) = -3aHc,/{2r), (40) 

valid for r <; nsco- At the radius r — nsco + Ar, we have itr-(Ar/risco)^ — —iaHcs/nsco- The sonic point {ur = —Cs) is 
at Ar ~ 2\/aHrisco, and Ur = — Cs/2 at Ar ~ 2\/2aHrisco- Thus itj,(rs) ~ Cs/I/e, with 

E 



iE=(^^J ~2V^^i7^=2\/a/3r», (41) 

where /3 = Cs/(rf2). This should be compared to the disc thickness H = /3r: depending on the value of a, both Le < H and 
Le > H are possible. 

The value and sign of Lc depend on the thermodynamical and radiative properties of the flow, and cannot be estimated 
in a simple way. It is reasonable to expect \Lc\ ~ rs. 

5.3 Reflectivity at the Sonic Point 

We can understand qualitatively the effect of the transonic boundary condition on the p-mode by calculating the reffectivity 
TZs of the inner boundary. 

Consider a density wave Sh oc exp(i k dr) in the wave zone rin = r^ < r < riL, traveling toward the inner disc boundary. 
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Figure 8. The wave reflectivity at the transonic point of the inner disc as a function of the parameter Ly,/H, for Vs / Lc = (the left 
panel) or Ts/Lc = —3 (the right panel). In both panels, the heavier lines are for Cs = O.lrf! and the lighter lines for Cs = 0.05rO. The 
short-dashed, solid and long-dashed lines are for m = 1,2,3, respectively. The wave frequency is set to be i^; = 0.7mQ(rs). 



Upon reflection, the wave becomes Sh cc TZs exp(— i kdr). Including the correct WKB amplitude (see Tsang & Lai 2008a), 
the wave outside the sonic point can be written as (up to a constant prefactor) 



Sh^A 



J rs 



exp til kdr \ + TZs exp ( —i kdr 



{rs <r < hl) 



where 



-D) 



1/2 



D 
rEfc 



1/2 



(42) 



(43) 



To apply the boundary condition (|36|) to equation (|42|) . we neglect Ur in equations (|6])-(|8]l at r = -I- e, with e <C Tj, and 
rs — risco- Implicit in this procedure is the assumption that the fluid perturbations do not vary signiflcantly between rs and 
rs + e. We then obtain 

ik + L^^ + K - 2mQ,/{ruj) 



Us 



where 



and 



ik-L-^^ - K + 2mQ./(r 



A' 



K = 



(44) 



(45) 



(46) 



l + {L^/r)-i{CoLT./2cs) 

When considering the damping of the p-mode due to the transonic flow, the quantity |7?.sP — 1 is the most relevant (see 
section 5.4 below). Let 

K = \K\eiLp{-iiP), (47) 

we have 

4,k\K\ simp 



\ns\-i = - 



{L^^ + \K\ cos 7/;)2 + {k + \K\ sin V')^ 



(48) 



^ Note that since the group velocity of the wave has opposite sign as the phase velocity for r < rjL, the wave of the form cxp(i kdr) 
(with A; > 0) is inward propagating. 
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Using /? = Cs/{rQ), Hi = — cimfiisco (where < < 1), we find from equation (146 p that 

V' = tanM^- SS-FT + tan M — -r , , , , ■ 49 

Figure 8 shows how the reflectivity depends on various parameters of the disc inner edge. In particular, for small Ls/H — 
L-E/{(3rs), i.e., when the surface density of the disc decreases rapidly at the sonic point, \TZs\^ is only slightly smaller than 
unity and the wave loss at the inner edge of the disc is small. 



5.4 Mode Growth Rate in the WKB Approximation 



Consider the p-mode trapped between rin = — rjsco and riL. With the reflectivity at riL given by TZ (see section 3), we 
can write the wave amplitude for r < riL afl 

1/2 r 



Sh oc 



D 
rT.k 



exp 



kdr ] +TZ exp I i 



k dr 



{ts < r < hl) 



(50) 



On the other hand, with the reflectivity at rin — given by TZs , the wave can also be expressed as (|42|l . For stationary waves 
we therefore require 



exp(2je) = 7^7^s, with e = /" kdr = er 
where Or and Oi are real. The real eigen- frequency cur is given by 

e,. : 



kr dr - 



dr — nil + 



Cs 2' 

where TZTZs = \TZTZs \ cxp{iip), and n is an integer. The mode growth rate cui is determined by |7?.7?.s| = exp(— 20i), or 

!7^7^J - r 



tanh Qi 

For Qi = J^'^ kidr <^ 1 and ki ~ ujiQr/icsV^'r ~ i^''), we obtain 

\\mzs\ + i 



7^7^J + 1 



(51) 



(52) 



(53) 



-dr 



(54) 

where we have assumed ubr <0. Equation (|54|l is to be compared with (|19p . where perfect reflection at rin is assumed. Clearly, 
to obtain growing modes we require \TZTZs\ > 1. For a given \TZ\ > 1, growing modes are possible only when the loss at the 
sonic point is sufficiently small (i.e., \TZa\ is sufficiently close to unity). 



5.5 Numerical Results 

We solve equations (|6])-(|8l) (with Ur = 0) subjected to the radiative outer boundary condition (|26p and the transonic inner 
boundary condition p6p . 

Figure 9 depicts an example of the p-mode wavefunctions. Again, the discontinuity in the angular momentum flux F at 
Tc signifies wave absorption; since rc < rpeak, this leads to mode growth. Comparing with Fig. 4, here the angular momentum 
fiux at Tin is significantly nonzero, indicating wave loss through the sonic point. Nevertheless, the corotational instability is 
sufficiently strong to overcome the loss and makes the mode grow. 

Figures 10-11 show the fundamental p-mode frequencies and growth rates as a function of the disc parameters. Consistent 
with the result of section 5.3 (see Fig. 8), growing modes are obtained for sufficiently small Ls- Negative Lc also tends to 
reduce wave loss at and make the growing modes possible. Such values of Le and Lc are not unreasonable for black hole 
accretion discs. 

It is important to note that while the mode growth rates uii depend sensitively on the inner disc parameters, particularly 
the physical property of the transonic flow near the ISCO, the real mode frequencies Ur show only weak dependence on the 
inner disc parameters (e.g., uir decreases with increasing sound speed; see Fig. 6). Thus we may expect that kHz QPOs appear 
only in certain accretion states of the black hole, and the frequencies do not vary much as the accretion rate changes. 



6 THE ROLE OF ROSSBY WAVE INSTABILITY 

Lovelace et al. (1999) (see also Li et al. 2000) have shown that when the vortensity — K^/(2f2E) has an extremum at a 
certain radius (rpcak) in the discjfl, it is possible to form normal Rossby modes around rpeak- If the trapped Rossby waves 

Note that this definition of 7?. diff'ers from that in Tsang & Lai (2008a) by a phase factor of cxp(j7r/4). 
^ Lovelace et al. considered non-barotropic flows, so the "generalized vortensity" depends on the entropy profile of the disc. 
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Figure 9. Wavefunctions for a disc p-mode. Tlie notations are the same as in Fig. 4. The disc has sound speed Cs = O.lrH and constant 
density profile (p = 0), and the m = 2 mode is calculated using the transonic inner boundary condition II36I I with L-^/H = 0.25 and 
Lc = oo. The eigenvalues are LOr = 0.725mSlisco a-nd uji/u)r = 0.00267. 



can propagate on both sides of the corotation, a standing pattern of waves of opposite energies are formed, making the 
mode unstable — This is the "Rossby wave instability" . Tagger & Varniere (2006) have considered the MHD version of the 
instability and suggested that it played a role in the diskoseismic modes around black holes (see also Tagger 2006). 

We do not find any trapped Rossby modes in our calculation. To clarify the issue in hght of works by Lovelace et al. and 
by Tagger & Varniere, let us consider equation (llOfl and define the effective potential 



, , 2mn / d , D 
Kff r = — :- — In — + — + _. 55 

rcj \dr D J r-' cj 

The wave equation can be approximated by {d^ /dr^ — Vcff)^^ — (see Tsang & Lai 2008a). We will focus on modes with rc 
very close to rpeak (i.e., \rc — r-pcaki < rc, so that ujr ^ mfipoak; see Fig. 1). For \r — rpcak] < r-pcak in a thin disc (so that 
rn^/r^ can be neglected compared to K^/cf), the effective potential becomes 

2mQ, f d \ 2 / r- rpeak \ , 

where in the second equality we have used Rc = rc — iircOJi/ quJr), ^ oc r"'', and defined Lc^ via 

dhlC ^ _r J^pcak ^ ^j,^^ _ rpcak] < rpeak). (57) 

Consider the case Vc < Tpcak and assume coi <^ ojr- The Rossby wave zone (where VeH < 0) hes between Vc and rc + Ar_R, with 

where H ~ Cs/k, and we have used L(^/H ~ r/ H 2> 1. The number of wavelengths in the Rossby zone is 

kdr^ (-Kff)'/' dr ~ 4^(''p-^k rc) ^^g^ 

Two points should be noted: (i) Since XT^+'^'^fedr « 1 for L( ~ rpeak, no stationary wave can form in the Rossby zone; (ii) 
since the Rossby zone lies only on one side of the corotation radius, even if the mode can be trapped it will not grow by the 
Rossby wave instability mechanism. Similar result can be obtained for the rc > rpeak case. We conclude that for the "smooth" 
vortensity maximum (with length scale ~ r; see the lower panel of Fig. 1) considered in this paper, there is no trapped 
Rossby mode around rpeak and the Rossby wave instability is ineffective. 

In the hypothetical situation where the vortensity has a minimum at r = rmin, equation (I56|l should be replaced by 



Kf.(0--^(^)+^, (60) 



© 0000 RAS, MNRAS 000, 000-000 



16 D. Lai and D. Tsang 




Figure 10. The real and imaginary frequencies of disc p-modes (with azimuthal wave numbers m = 1, 2, 3) as a function of Ly,/H [see 
eq. Il35l l] The modes arc calculated using the transonic inner boundary condition I I36I I with Lc = oo. The disc has a constant surface 
density profile and the sound speed is Cs = O.lrQ, (left panels) or O.OSrn (right panels). 




Figure 11. The real and imaginary frequencies of disc p-modes (with azimuthal wave numbers m = 1,2,3). The modes are calculated 
using the transonic inner boundary condition (1361 . The disc has a constant surface density profile and the sound speed is Cs = O.lrC 
The left panels show the cases with = —rg/3 by varying L-^/H, and right panels show the cases with L-^/H = 0.25 by varying Ts/L^ 
[see eq. | |35| |]. 
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where we have used 

c? ln(^ r — Vi, 



dr Ll ' 



(for |r 'min 'min ). (61) 



In this case, for a mode with ujr = mn(rmin) (or Vc = rmin), we find Veff(rc) ^ — 2/(gL^) + (for a»i <C Wr-)- When 

Ve« {re) < 0, or when 

Lc < - H, (62) 



Rossby waves can propagate on both sides of the corotation, leading to mode growth this is the Rossby wave instabihty. 
Thus, the Rossby wave instabihty would operate if there existed a "sharp" vortensity minimum in the disc (with ^ varying 
on the lengthscale comparable or less than the disc thickness) — this is not the case for typical black hole accretion discs 
considered in this paper. 



7 DISCUSSION 

High-frequency QPOs (HFQPOs) in black-hole X-ray binaries have been studied observationally for more than a decade now 
and they provide a potentially important tool for studying the strong gravitational fields of black holes (see Remillard & 
McClintock 2006). Despite much theoretical effort, the physical mechanisms that generate these QPOs remain unclear (see 
section 1.1 for a brief review of existing theoretical models). Ultimately, numerical simulations of realistic accretion discs 
around black holes may provide the answer. However, such simulations are still at their early stage of development and have 
their own Umitations (e.g., De Villiers & Hawley 2003; Machida & Matsumoto 2003, 2008; Arras et al. 2006; Fragile et al. 2007; 
Reynolds & Miller 2008; Beckwith et al. 2008; Shafee et al. 2008; Noble et al. 2008), semi-analytical study remains a useful, 
complementary approach in order to identify the key physics involved. 

In this paper, we have studied the global instability of the non-axisymmetric p-modes in black-hole accretion discs. These 
modes have frequencies to ~ (0.5 — 0.7)mSlisco (where m is the azimuthal wave number, Slisco is the disc rotation frequency 
at the inner-most stable circular orbit), where the pre-factor (0.5-0.7) depends on the inner disc structure. Recent works (Arras 
et al. 2006; Reynolds & Miller 2008; Pu & Lai 2008) suggested that, unlike other diskoseismic modes (g-modes and c-modes), 
the p-modes may be robust in the presence of disc magnetic fields and turbulence. Our linear analysis showed that due to 
GR effects, the p-modes may grow in amplitude due to wave absorptions at the corotation resonance. For a given m, only 
the lowest-order p-mode has sufficiently high frequency {ui > mOpeak; see Fig. 1) to be driven overstable by the corotational 
instability, while high-order (lower frequency) modes are damped by the corotational wave absorption. 

The greatest uncertainty of our calculation of the p-mode growth rate concerns the boundary condition at the inner disc 
edge near the ISCO. In particular, the rapid radial inflow at the ISCO has the tendency to damp the mode (see Blaes 1986). 
While our analysis in section 5 indicates that this damping does not completely suppress the mode growth under certain disc 
conditions, it suggests that mode growth may not always be achieved in real black-hole accretion discs. Observationally, it is 
of interest to note that HFQPOs are observed only when the X-ray binaries are in the steep power-law state, while they do not 
appear in other spectral states (Remillard & McClintock 2006). In particular, HFQPOs are absent in the thermal (soft-high) 
state, believed to correspond to geometrically thin discs extending down to the ISCO. It is reasonable to expect that in this 
state p-modes are damped due to the rapid radial inflow. 

Our current understanding of the steep power-law state (also called very high state) of black-hole X-ray binaries is rather 
limited. A thermal-radiation-emitting disc is suggested by spectral modelings, but it is not clear whether the disc is truncated 
at the ISCO or slightly larger radius (see Done et al. 2007). The observed power- law radiation component requires a signiflcant 
corona that Compton up-scatters the disc thermal radiation. It is possible that in the steep power-law state, the inner disc 
behaves as a more reflective boundary (modeled in section 4) than a transonic flow (modeled in section 5), and thus more 
robust p-mode growth can be achieved. One possibility is that a significant magnetic field flux can accumulate in the inner 
disc when the disc accretion rate is sufiiciently high (see Bisnovtyi-Kogan & Lovelace 2007; Rothstein & Lovelace 2008 and 
references therein). Such a magnetic field may also enhance the corotational instability and induce variability in the power-law 
radiation fiux (see Tagger & Varniere 2006). 

Although the p-mode growth rates depend sensitively on a number of (uncertain) disc parameters (particularly those 
related to the inner disc boundary), the mode frequencies axe more robust (see Figs. 5-7, 10-11). More precisely, the real mode 
frequency can be written as cur — CjmQisco, where a; < 1 depends weakly on m and has only modest dependence on disc 
parameters (e.g. sound speed). This implies a commensurate frequency ratio as observed in HFQPOs (note that in some of 
our models, the m = 2, 3 modes have the largest growth rates; see Fig. 5). The fact that a) < 1 would also make the numerical 
values of the p-mode frequencies more compatible with the measurements of the QPO frequencies and black hole masses. 

We note that our calculations in this work are done with a pseudo-Newtonian potential. For direct comparison with 



© 0000 RAS, MNRAS 000, 000-000 



18 D. Lai and D. Tsang 

observations a fully general relativistic calculatior[f] is needed including a careful treatment of the corotation singularity. 
Including the effect of black hole spin would likely increase the value of uj by modifying risco and fiisco, while uj will likely 
remain similar to the non-spinning case discussed above. We plan to study these effects in future work. 
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